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Abstract 

We propose a scheme for controlling interactions between Rydberg-excited neutral 
atoms in order to perform a fast high-fidelity quantum gate. Unlike dipole-blockade 
mechanisms already found in the literature, we drive resonantly the atoms with a 
state-dependent excitation to Rydberg levels, and we exploit the resulting dipole- 
dipole interaction to induce a controlled atomic motion in the trap, in a similar way 
as discussed in recent ion-trap quantum computing proposals. This leads atoms to 
gain the required gate phase, which turns out to be a combination of a dynamic 
and a geometrical contribution. The fidelity of this scheme is studied including 
small anharmonicity and temperature effects, with promising results for reasonably 
achievable experimental parameters. 
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1 Introduction 



The quest for the reliable implementation of a fundamental quantum gate is 
one of the central topics in the current research in quantum computation. 
The basic requirement is the controlled coherent dynamics [1] of the system 
realizing the quantum bits, while avoiding decoherence due to fluctuations of 
the external control parameters and coupling to an environment 1 . 

1 We dedicate this publication to Bruce Shore as one of the pioneers in developing 
the theory of coherent interactions between atoms and laser light. 
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In the last years, several proposals for different physical systems have ap- 
peared, ranging from trapped ions [2], photons in cavity QED [3], molecules 
[4], to quantum dots and Josephson Junctions [5]. A particularly attractive 
perspective is the possibility of using cold neutral atoms [6], whose steadily 
improving experimental control looks quite promising. Regarding the required 
logical operations, different possible universal sets of quantum gates have been 
devised, showing that single-qubit and two-qubit gates are sufficient to realize 
any iV-qubit unitary operator. In this context, a widely studied fundamental 
gate is the two-qubit controlled-phase gate, whose truth table is given by 
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where <f> is the phase and \g), |e) are the logical states. When equals n, 
this gate is equivalent to a Controlled-NOT and is universal if combined with 
arbitrary single-qubit rotations. 

An important figure of merit to evaluate a specific physical implementation is 
the ratio between the gate operation time and the coherence time characteristic 
of the system. In this sense, for a given quantum memory coherence time and 
gate error rate, faster gates are of course highly preferable. However, gate speed 
is limited by several system-specific factors, for instance by the characteristic 
energy scales that are present: roughly speaking, if the phase in Eq. (1) is 
produced by a state-dependent energy shift AE, the time it takes for a phase 
7r to be accumulated has to be of the order of AE' 1 . 

In trapped atomic systems the relevant energy scales are the frequency of 
the trapping itself (which can hardly exceed, say, a few hundred kHz) and 
the energy of atom-atom interactions. This has led to several theoretical ex- 
plorations towards employing stronger interactions like dipole-dipole forces 
between Rydberg-excited atoms [7] or molecular interactions giving rise to 
Feshbach resonances [8]. The main issue to be dealt with in this case is that 
strong interactions tend to perturb significantly the trapped-atom dynamics, 
thereby representing a potential source of errors. Ways around this problem 
depend again on the specific feature of the system under consideration. For 
instance, putting Rydberg atoms in a sufficiently stable and intense static 
electric field, a dipole-blockade mechanism can arise that shifts the two-atom 
interacting state out of resonance from single-atom exciting lasers [9,10]. In 
this two-qubit gate can be performed without ever switching on the 

dipole-dipole interaction, thereby avoiding decoherence due to entanglement 
of the qubit with the motion. 



2 



An alternative way to overcome trap-limited gate time scales is to employ 
quantum control techniques, as originally proposed in the context of ion-trap 
quantum gates in [11]. In this approach, the gate phase is not directly or 
indirectly obtained from a state-selective energy shift due to interaction, but 
rather by driving the two-particle system on a closed path in phase space, 
whose final result is to bring it back to its original state, while at the same 
time imprinting on it the desired phase. 

The present paper represents an application of these methods, developed in 
[11,12] for trapped ions, to the case of Rydberg-excited atoms already dis- 
cussed in [7]. The use of quantum control techniques allows for exploiting the 
strong perturbation due to dipole-dipole interaction as a means for driving 
the system into the goal state with a high fidelity that does not depend, to 
a good approximation, on the initial motional state of both atoms inside the 
trap, and in a time that is a fraction of the trap period. Moreover, since the 
phase is not accumulated as a result of some interaction-induced energy shift 
as discussed above, the atoms are not required to remain in the interacting 
state for the whole duration of the gate, but just during much shorter times 
as needed for imparting the "kicks" which drive their motion. In this way 
the effect of decoherence due to spontaneous emission from the highly excited 
Rydberg states is further reduced. 

A generic quantum computation requires in principle the ability to perform 
more complex sequences of two-qubit gates between any pair of qubits. To 
achieve this, one will need to move atoms between different locations in the 
quantum register. Thus the overall speed of a complete quantum computation 
in general will be basically limited by the time needed to transport the atoms 
close to one another before performing each actual gate operation on the de- 
sired qubit pairs. This limitation affects all quantum computing implementa- 
tions based on atomic (or ionic) systems. However, it can be partially overcome 
by applying quantum control techniques to the transport process [13]. In ad- 
dition, individual addressing, which is still an issue in present experiments, is 
not needed at the level of a two-qubit gate within the protocol proposed here, 
and therefore atoms can be kept during the gate at a shorter distance than 
required in other schemes by single-atom laser diffraction-limited addressabil- 
ity. Thus atom transport would only be required for two-qubit gates involving 
non-neighboring qubits, and it may be employed as well whenever a gate on 
a single qubit is needed during the computation, to single it out by moving it 
away from the other qubits. 

The structure of the paper is as follows. In Sec. II we introduce in general 
terms the simple control procedure we plan to apply to the system for this 
purpose; in Sec. Ill we describe the physical model for the implementation of 
the control pulses with dipole-interacting atomic systems and discuss exper- 
imental parameters; and in Sec. IV we analyze quantitatively the obtainable 
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Fig. 1. Phase space trajectory of a simple kick protocol: a) usual reference frame; 
b) rotating reference frame. The harmonic oscillator ground state is first given a 
positive momentum p\ = y/2px; then it is left evolve in the pure harmonic oscillator 
potential for a time 9 (in harmonic oscillator units), and hence kicked again by a 
negative momentum p2 = v2p2 = ~~ 2pi cos(#). The purely harmonic evolution time 
9 corresponds to the angle shown in the figure. The final kick, which brings the 
system back to the initial coherent state, has the same intensity pi as the first one. 

fidelity, before drawing our conclusions in Sec. V. 



2 Control scheme 



We begin by discussing a possible 'kick' scheme to manipulate the system 
phase while bringing the atoms back to the initial motional state, in the spirit 
of Ref. [11]. Although in the next sections we will adapt this protocol to our 
particular configuration, we point out that this scheme is applicable to quite 
general systems. In this section, the 'system' is simply given by a single particle 
in a one-dimensional harmonic potential, corresponding in the following to the 
relative motion of the two atoms. 

The basic idea is to lead the system along a closed path in the oscillator 
phase space by taking advantage of the momentum shifts - the so-called kicks 
- induced by some external force. The simplest possible protocol, i.e., that 
requiring the minimal number of kicks, is illustrated in Fig. 1 for the special 
case of the oscillator ground state. This scheme, as far as the kicks are ideal 
shifts (see the impulse approximation discussed in the next section), has the 
remarkable advantage of being independent of the initial state. 

We consider the unitary phase space displacement operator D(a) = e aat - a * a ; 
where and d are respectively the creation and annihilation operators of the 
harmonic oscillator and dimensionless units are used 2 . With this formalism, 
a kick of momentum p is simply given by e ipx = D(ip), where p = p/ a/2 is a 

2 We recall that D(a)\0) = \a), where |0) is the oscillator ground state and the 
coherent state \a) is an eigenstate of the annihilation operator a\a) = a\a). 
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real number and x = (a + a))/ \/2 is the position operator. The protocol shown 
in Fig. 1 is a sequence of three kicks alternated by free harmonic evolution, 
the latter being governed by the harmonic time evolution operator Uho(t) = 
e -i(at&+i/2)t_ By using the well k nown relations D{a)D{(3) = D(a + P)e iIm ^ 

and Uh (t)D(a) = D(e~ lt a)UY lo (t) ) one can easily rewrite the kick operator 
G = D(ip 3 )U ho (02)D(-ip 2 )U ho (di)D(m) in the form 

G = U ho (0i + ^ 2 )£>(i[pi - P2e Wl + p^+^e^ , (2) 

where gcO m = P1P3 sin(#i + 6 2 ) — pip 2 sin 6>i — p 2 pz sin 6 2 . The 'triangle' scheme 
is nothing but the condition p 1 -p 2 e Wl + p 3 e i(ei+f?2) = 0, so that G = U ho (6 l + 
2 ) e i0 g com an( j ^he initial motional state is restored. The triangular path is 
evident in the rotating frame description (see Fig. 1(b)), namely the interaction 
picture where the action of the oscillator Hamiltonian = a) a +1/2 is 
made implicit and the gate operator does not contain the term £/ho(#i + #2)- 
In addition, one immediately finds pi^sin^ = p 2 P3sin# 2 = PiPz^-iPi + 
6 2 ) = 2 A, where A is the area of the triangle described by the rotating phase 
space trajectory. Finally, gcO m — —2A and, for the considered case < 61 = 
9 2 = 9 < 7r/2, 0g CO m = — Pi sin(26 ) ). Being based on operatorial identities this 
result is evidently independent of the initial state. It is worth recalling the 
'geometric' nature of such phase, as discussed in Ref. [12], reflected by the 
non-commutativity of the phase-space displacement operators. 

We stress that the discussion in this section simply refers to the dynamics 
of a ID harmonic oscillator. Below we will see that such scheme can be ap- 
plied to the effective ID oscillator obtained for the relative motion of the two 
considered atoms, so that this geometric phase will be used in a protocol for 
a two-qubit gate. The concrete implementation will however introduce some 
approximations whose validity will actually depend on the initial oscillator 
state (see Subsec. 3.3). 



3 Physical implementation 



We now turn to the analysis of the particular case where the kicks are provided 
by the switching of the dipole-dipole interaction. We will first give an overview 
of the actual form of the interaction, then we will discuss its application to 
the protocol presented in the previous section, and finally we will examine the 
validity conditions for the approximations introduced. 
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Electric field 

Fig. 2. Level scheme for the quantum phase gate. The internal state \g) is coupled 
by the Rabi frequency to the Rydberg state |r). Since the state |e) is unaffected 
by the laser, the dipole interaction arises only if both atoms are initially in the \g) 
state. The harmonic trapping potential Vj r ) felt by the atoms in the Rydberg level 
|r) is different from the potential V\ g \ i e \ of the hyperfine levels |e) and \g). 

3.1 The dipole- dipole interaction 

We consider two atoms in two different one- dimensional harmonic wells with 
the same frequency u, separated by a distance xq. This system can be thought 
of as an approximate description of two atoms sitting in different sites of a 
optical lattice or in distinct trap minima on an atom chip, where the dynamics 
transverse to the direction joining the atoms are completely frozen. We will 
come back on such an approximation in Subsec. 3.3. Each atom corresponds 
to one qubit, where two internal hyperfine atomic levels, which we call \g) and 
|e), represent the logical states. 

Assuming an underlying static electric field, the gate is performed by driving 
state-dependent transitions to Rydberg levels (see Fig. 2), switching on and off 
a dipole atom-atom interaction which allows to manipulate the global phase 
of the system, as described in Ref. [7]. Initially, each atom can be either in 
\g) or |e). The internal state \g) is coupled by a laser to a Rydberg state |r), 
where the transition \g) — > \r) is assumed to be obtained via a Doppler-free 
two-photon absorption, which suppresses any photon kick contribution and 
gives rise to an effective Rabi frequency fl (see Fig. 2). The state |e) is instead 
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fully decoupled from the Rydberg manifold. In this way, we switch on a dipole- 
dipole interaction only when both atoms are initially in the internal state \g). 
This interaction is kept on for a very short time, after which the atoms are 
brought back to the \g) states. While in the \rr) state, the atoms acquire a 
dynamical phase as well as a momentum shift (kick). Thus, repeating this 
procedure according to the scheme described in Sec. II, such a process realizes 
the phase gate (1), provided we undo the 7r-phase due to the Rabi flopping for 
the states \ge) and \eg). 

The Hamiltonian of the system without the laser coupling between the states 
\g) and |r) reads 



i=l,2 L 



E 



=1,2 



Pi 2 / \ 2 

1 — muj (Xi — xoi) 

2m 2 v ; 

Pi /2 / \1 

h -mw [Xi - xoi) 

2m 2 



(\9)M + \e)M) 
\r)i(r\ 



(3) 



+ V dip {x 1 - x 2 )\r) 1 (r\ <g> \r) 2 (r\ 



where we consider that the (harmonic) trapping potential is the same for 
the two hyperfine levels (frequency uj) and different for the Rydberg state 
(frequency uj'). The general 3D expression for the dipole-dipole interaction is 



V dip (r) 



1 



4vre 



(4) 



where fti, fi 2 are the dipole moments of the two atoms, which depend on the 
Rydberg level and the external electric field, and the 3D relative distance r 
reduces in our case to the ID relative motion coordinate x — x 2 — x\. In the 
Hamiltonian (3) m is the mass of the atoms and we take x 02 — x i = x , where 
x is the distance separating the two harmonic wells (see Fig. 2). 

The possibility of exploiting this interaction to gain a non-trivial phase shift 
and effect a quantum gate has been suggested in Model A of Ref . [7] . In this 
scenario, during the interaction time At, the system acquires a dynamic phase 
uAt, where u is the energy shift created by the dipole potential. However, in 
order to reproduce the truth table (1), one also has to ensure the restoring of 
the initial motional states, which are strongly affected by the non-negligible 
mechanical effects caused by the dipole interaction - and this is not straight- 
forward in this scheme. To overcome this, in Ref. [7] a Model B has been 
introduced, where the strong dipole-dipole interaction is instead used to in- 
duce a dipole-blockade mechanism that still ensures state-dependent dynamics 
while actually avoiding mechanical back-action on the atomic motion. 



An alternative solution can be achieved by controlling such strong mechanical 
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effects to gain an additional geometrical phase, while at the same time bringing 
the system back to the initial motional configuration, similarly to the method 
proposed in Ref. [11]. 

In the following we will consider only two possible configurations for Vdi P (x): 
the repulsive case V + (x) = /ii/X2/47re a; 3 , corresponding to dipoles aligned 
in the direction orthogonal to the and the attractive case V-(x) = 

— 2fj J ifj J2 /4:7ie x 3 = — 2V + (x), corresponding to dipoles parallel to the x-axis. 
This in principle could be achieved by rapidly changing the electric field ori- 
entation between subsequent laser pulses - although from a practical point of 
view it may be problematic to realize this on the time scales required here. A 
viable alternative would be to keep the external field fixed and excite the sys- 
tem to Rydberg states having differently oriented dipole moments, something 
that could be performed without requiring individual laser addressing of the 
atoms. 

The expression of the dipole moment for Rydberg states is given by = 
(3/2)ni i 2q , i,2eao [14], where e is the electron charge, a is the Bohr radius, 
and n, q are quantum numbers which label the different states. The quantum 
number q obeys the relation q = n — 1 — \m\,n — 3 — \m\, —{n — 1 — \m\), 
where m is the magnetic quantum number. Typical values considered here (for 
87 Rb) are m )2 = qi,2 + 1 = 50. 

We prefer to work with dimensionless quantities by using the oscillator units 
of the harmonic potential with frequency u>, felt by the atoms for most of 
the time. The dipole potential in particular can be rewritten as V + /hcu = 
a+/(^/aho) 3 , where a\ = (9/4)nigin2g2(/V m e)/(o4io/ao), A* = m /2 is the re- 



oscillator length and then we put H — [i = a\ iQ = 1. 

Note that in order to realize an ideal transition to the Rydberg state one needs 
the condition Q » \V±(xq)\, i.e., Q ^> 2a 2 + /x^, which is easily satisfied in all 
the cases of practical interest. Realistic values for the parameters are discussed 
in Subsec. 3.3. 

We describe the action of the dipole potential arising when both atoms are in a 
Rydberg state by using an impulse approximation. We indeed assume that the 
dipole interaction is very strong compared to the harmonic trapping and that 
its duration is very short, so that the system undergoes a small displacement 
during this time and the corresponding harmonic potential contribution to the 
phase is negligible. Furthermore we separate the centre of mass and relative 
motion, where only the latter is affected by the dipole-dipole interaction. Due 
to the very short time spent by the atoms in the Rydberg state, we then 
consider the centre of mass motion to be completely unperturbed by the gate 
and will only examine the relative motion. For a short time interval At the 




harmonic 
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relative motion time evolution operator can hence be approximated by 

U At ~ e - iV ^ x)At , (5) 

where, in force of the impulse approximation, we dropped both the kinetic 
term and the harmonic oscillator potential. In addition, we expand the dipole 
potential around the centre Xq of the harmonic trap 

V dip {x) ~ Vdip(xo) + V^ ip (x )(x - x ) , (6) 

where the prime denotes the first spatial derivative. The largest possible value 
of \x — xo\ is fixed by the extension of the relative motion wavefunction ip. This 
can be quantified by the average position (x) = J x\ip\ 2 dx and the mean square 
length (Ax) = ^ (x 2 ) — (x) 2 , so that the validity of the linear approximation 
(6) requires 3 \x — x \ < \(x) — x \ + (Ax) <C x . A full description of the 
consequences of the above conditions in terms of experimental parameters is 
given in Subsec. 3.3. 

With the previous approximations the time evolution operator acquires the 
elementary form 

U At ~ e - iuAt + iA p(x~*o) ? (7) 

where u = Vdi P (a;o) is the zero-energy shift and the momentum transfer ('kick') 
Ap is simply given by the classical expression Ap = — V'(x )At. 



3.2 The dynamical phase 



The 'ideal' picture described in Section 2 can now be joined with the actual 
form of the kick dynamics. The phase space depicted in Fig. 1 will clearly 
correspond to the relative motion harmonic oscillator with frequency uo (which 
scales to 1 in dimensionless units), whereas the kicks will be provided by the 
dipole potential. As previously discussed, in first approximation one has p\ = 
\V+(xo)Ati\ and P2 = |V^.(xo)At 2 |, i.e., the momentum transfers p± t 2 require 
a finite (although very short, as far as Vdi P is large) time Ati t2 - In addition, 
since p 2 = 2p 1 cos# and 4 VL(xo) = — 2V r _[(xo), one has At 2 = Aticos^. 
The total phase gained by the system is the sum of the dynamical phase 

3 One could improve this approximation by linearizing Vdi p (x) around (x) instead 
of xq. This, however, would make the calculation dependent on the initial motional 
state. A proper discussion of the corresponding procedure could be done within the 
optimal control technique. 

4 Here we assume that all the kicks are obtained by coupling the atoms to the same 
Rydberg state, choosing the simplest experimental situation where only one laser 
frequency is required. In principle, nothing prevents from using more complicated 
schemes which add further degrees of freedom to the protocol. 
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0d y n = -[V + (x )At 1 + V-(xo)At 2 + V + (x )At 1 \ = -2V + (x )At 1 (l-cos9) and 
the geometrical phase geO m- Hence 

a 2 Q rv 4 

= -2^Afi(l-cos0)-^A*?siii(20) , (8) 

which can be solved to give 



a\At x = 



(9) 



Xr, 



9 sin(2#) 



cos 9 — 1 



1 



(cos#- 1) 



9 sin(2#) 



Xq 



This allows to choose the proper combination of Rydberg quantum numbers 
and interaction time in order to achieve the desired phase. The full gate time 
required by the operation will be r = Aii(2 + cos#) + 29 ~ 29. It is also 
interesting to consider the fast quantum gate limit 9 — > where one easily finds 

-(x /3)9 3 /\ 



a 



Ati ~ (xq/3)J\(P\/9. This implies 0dyn 



which means that in this limit the global phase 
geometrical contribution. 



and 



y gcom 



is completely given by the 



Let us stress again that the phase <j) will arise only if the atoms are initially in 
the \gg) internal state, since otherwise the dipole potential is never switched 
on. 



3.3 Validity of the employed approximations 



Our impulse approximation (5) requires that, during the interaction time At, 
the action of the relative motion oscillator Hamiltonian H = p 2 /2 + (x— x ) 2 /2 
(we keep on using harmonic oscillator units) be negligible 5 . The contributions 
of the different terms in the Hamiltonian can be quantified at a given time by 
taking their expectation values on the corresponding state \ip{t)). We therefore 
require 

W<I0W)| (io) 

during the whole kick time At. However, from the point of view of the gate 
dynamics, it is also necessary to have 

{H ) At « |0| , (11) 

since otherwise we should include such a contribution in the calculation of the 
dynamical phase 0d yn - 

5 In this analysis we assume for simplicity that u and u' are of the same order [7] . 
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Besides conditions (10) and (11), we also have to consider the validity require- 
ments of the linear approximation (6), i.e., the condition 



(x) - x \ + (Ax) < x , 



(12) 



necessary for the linearization of the dipole potential. Clearly, the explicit form 
of all the previous conditions depends on the states reached by the system 
during the dynamics of the kick scheme. 

For example, if the system stays very close to the oscillator ground state 
during the whole gate time, one can estimate the above expectation values 
by using \ip{t)) = |0). In this case (H ) = 1/2, (x) = x and (Ax) = l/y/2. 
Linearization then requires x ^> l/y/2, while, by assuming the validity of the 
latter condition, Eq. (10) reduces to 1/2 <C V + (x ), which can be made explicit 
by evaluating V + (x ) = a\/xl with Eq. (9). Eq. (11) is simply At < 2|0|. 

Since we want our protocol to be independent of the initial motional state, 
it is interesting to find the conditions under which our approximations are 
applicable to a generic state not too far from the ground state. We therefore 
consider the case \tp(t)) = \cx(t)), where \a(t)) is a coherent state subject to 
the constraint \cx(t)\ < R, R > 0. Then, if the initial state of the system is 
the coherent state |ck(0)), the boundary R can be estimated within the ideal 
kick scheme, where \cx(t)\ < |«(0)| +p\. Concerning the expectation values we 
have (H ) = \a\ 2 + 1/2 < R 2 + 1/2, \(x - x )\ = v^lRea] < V2R, and still 
(Ax) = 1/ y/2. If Xo ^> y/2R +1/ V2 then the use of the linear approximation 
is legitimate and Eq. (10) becomes R 2 + 1/2 <^ V+(xo), while Eq. (11) is 
(i? 2 + l/2)At<2|0|. 

We explicitly evaluate these conditions for the most interesting case a(0) = 0, 
where initially the system is in the oscillator ground state and hence R — pi — 
(3/\/2)V + (xo)Ati/xo- Qualitatively, fast gate performances require small val- 
ues of 9, while the linear approximation (6) needs large values of xo- However 
large distances disfavour the impulse approximation (5), which is valid for 
strong dipole interactions. High Rydberg quantum numbers are then neces- 
sary to enhance the dipole moments. In order to study the interplay among 
the different parameters it is useful to consider Eq. (9) in the limits 9 — > and 
xo — > oo. Then, it turns out that in all the cases of practical interest condition 
(10) is already implied by Eqs. (11) and (12), so that one is left with only two 
conditions. In addition, the most stringent form of these requirements takes 



place in the 9 — > limit, where V + (x )Ati ~ (x /3)y \4>\/9. In conclusion, the 
linearization condition (12) and Eq. (11) respectively reduce to 





(13) 
(14) 



At! < 29 , 
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valid in the fast-gate limit 6 — > 0. 

The conditions (13-14) are actually rather demanding in terms of experimen- 
tal parameters, even though certainly achievable. Assuming Eqs. (13) and 
(14) to be well satisfied by a factor of ten difference between the two terms 

of each inequality so that x = lO^/|</>|/0 and Ati = 9/5, and recalling the 

9 -> result a\ ~ (arg/3A*i)^/|0|/0, one gets a* ~ (lO 5 /6#)(|0|/#) 5 / 2 . Un- 
fortunately a+ cannot be arbitrarily large. For example, for the case of 87 Rb 
atoms, where fx = 87/2 amu, with a harmonic frequency uo = 2n x 100 KHz, 
so that (fx/m e )/(aho/cLo) = 86.8, by choosing the Rydberg quantum numbers 
ft<i,2 = 91,2 + 1 = 50 one has a+ = 1.2 x 10 9 . Then, for |0| = n (equivalent to 
the C-NOT gate), the shortest gate-time obtainable within this rough upper 
bound for a\ corresponds to 9 = [{10^p\) 5 /6a 2 + ] 2 / 7 = 0.09, i.e., r ~ 0.03T ho , 
where Th = 2ir is the harmonic oscillator period. Notice also that, while 
the condition corresponding to Eq. (14), Ati = 9/5, can be easily satisfied 
experimentally, Eq. (13) for such a small value of 9 imposes a very large sep- 
aration between the harmonic wells (in harmonic oscillator units), namely 
x = 10^J\4>\/9 = 58, significantly larger than the value corresponding to the 
typical distance between neighbouring sites in currently employed optical lat- 
tices (for example, for a laser wavelength of 800 nm, the above parameters 
yield xq = 8.3). However, this requirement can be clearly relaxed by increas- 
ing the gate time (see the next section). Besides, there is the possibility of 
increasing xo independently of the lattice spacing, e.g., by using superlattices 
or the methods discussed in Ref. [15]. 

In the final part of this section, we briefly discuss the validity of the one- 
dimensional approximation (see also [7]). We assume a harmonic trapping 
also in the radial direction, with a frequency uj± and restore the dimensional 
units. The radial degrees of freedom could then be neglected if the level spacing 
fkv± were much larger than all the other energy scales which come into play. 
Clearly, this is difficult to accomplish as far as one also needs large values of the 
frequency uo along the x direction, necessary to increase the ratio £o/ a ho y/ui. 
However, provided each atom is sitting in the ground state of the corresponding 
harmonic well, the dimensionality does not play any role. Indeed, assuming 
the exact dipole polarization, the momentum transfer takes place only in the 
x direction and the gate dynamics is decoupled from the radial degrees of 
freedom, irrespectively of the value of u±. In the case where the initial motional 
state significantly differs from the ground state, one instead has to compare 
hcu± with hio and with the component of the dipole potential affecting the 
radial motion, which is however smaller and smaller as the inter-well distance 
increases. Consequently, the validity condition of the linear approximation also 
favours the quenching of 3D effects, even though a quantitative analysis of this 
approximation would require a full description of the 3D dynamics, which is 
beyond the scope of this paper. 
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4 Fidelity analysis 



We are now interested in considering possible error sources. The effects we take 
into account in the following are (i) contributions not included in the impulse 
and kick approximations, which are important when the validity conditions 
discussed in the previous section are not well satisfied, and (ii) anharmonic 
terms in the trapping potential. Notice that, since as far as our approxima- 
tions are valid the model is independent of the initial motional state, finite 
'temperature' (in a sense specified later) effects are actually included in the 
first point. 

Note that in addition, as discussed in [7], minor error sources that are com- 
pletely neglected here are the spontaneous emission and the black body radi- 
ation for the Rydberg state. 

In our scheme atoms in different wells have a negligible overlap and can be 
considered distinguishable. In this case the gate fidelity can be defined as [16] 

F = min{Tr ext [(x|W T (| X )(x| ® Po)Z4|x>]} , (15) 

where U T is the unitary operator governing the evolution of the full system 
during a single gate operation, \x) = J2 a /3 c apWP) is an arbitrary internal state 
of the two atoms with the normalization J2 a /3 \ c a/3\ 2 = 1, |x) is the state given 
by the application of the ideal gate (1) to \x), and po is the initial density 
operator relative to the external degrees of freedom, over which the trace is 
taken. In our case, the external states are simply given by the motional states. 
Furthermore, since the centre of mass motion is unaffected by the gate, we 
can restrict the analysis to the relative motion 6 . 



4-1 Fidelity for pure motional states 

Let us consider the particular case p — \k)(k\, where \k) are the eigenstates 
of the relative motion harmonic oscillator. This corresponds to assuming zero 
entropy, the system being in a pure state. 

We can write the action of U T as 



6 This is not generally true, e.g., if the atoms were in the internal state \er) for a 
time comparable with the oscillator time scales. 
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Fig. 3. Upper panel (a): fidelity F = (l + Re(^id|V'num))/2 as a function of the time 6 
and the distance xq (harmonic oscillator units), for \(f>\ = ir and At\ = 10 -4 x2ir. The 
state iV'num) is obtained by calculating numerically the full gate-dynamics (in the 
absence of anharmonic perturbations) for the state \gg) ® |0), while the ideal state is 
simply \1p1A) = e*^|0). The contour lines drawn in the lower plane correspond to the 
fidelity values 0.9,0.99,0.999, while the thick line plots the relation x® = W^/ir/9. 
Lower panel (b): same as (a) but in the presence of an anharmonic quartic term 
with A = 10~ 2 . 

U T {\X) ® \k))=Ylc aP \a(3) ® \k)e^ k+1 ^ + (16) 

a/3 

+c 99 \gg)®Y. a ^\k')e-^' +l l^ , 

k< 

where the primed summation does not contain the term af3 = gg, and J2w \ a kk'\ 
1. The explicit calculation then gives 



F k = min{l + 2(1 - \c gg \ 2 )\c gg \ 2 [Re(e^al k ) - 1]} = 

\ c 99l 

_ 1 + \a kk \ cos(4> - <p kk ) 
2 

where 4> kk = &rga kk and e l ^a* kk = (V ; num(^)|V ; id(^)) is the overlap between 
IV'id(fc)) — e^\k) and the final motional state obtained by applying the full 
gate dynamics to the state \gg) ®\k). 

We compute the zero-entropy fidelity in the good initialization case k = 0. In 
the ideal case where the impulse approximation is valid one has «oo — e 1 ^ 00 = 
e 1( ^ and a^ k — for k > 1. 

In Fig. 3(a) we show the gate fidelity F when the protocol is applied to the 
motional ground state (fc = 0), as a function of the two most important pa- 
rameters 8, Xq. This allows a direct comparison with the conditions discussed 
in Subsec. 3.3. The data are obtained by solving numerically the time depen- 
dent Schrodinger equation corresponding to the full kick dynamics applied 
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to the oscillator ground state, for the case of purely harmonic trapping. The 
behaviour corresponding to Eq. (13) is also plotted, showing a qualitative 
agreement with the curves at constant fidelity. Indeed, for the considered pa- 
rameters, Eq. (14) is always very well satisfied and the terms not included in 
the linear approximation are the main reason of discrepancy with respect to 
the ideal scheme. 

In Fig. 3(b) we repeat the same analysis by adding to Hq an anharmonic 
term — \{x — x ) 4 /2, corresponding to the first correction one would have for 
the square sinus potential of an optical lattice 7 . Only a very slight worsening 
of the fidelity is observed for the realistic value A = 10~ 2 , showing that in 
practice the validity of the linear approximation is the crucial condition. 



4-2 Finite temperature and entropy effects 

In view of a realistic quantum computing scenario it is important to investigate 
the dependence of the gate performances on temperature. Indeed, in the course 
of the computational process the system is expected to undergo a certain 
heating - for example due to moving the atomic qubits around by shifting the 
optical lattice - which could in principle be limited by a periodic re-cooling of 
the qubits [17]. 

On the other side, atoms in optical lattices can be considered practically iso- 
lated from the environment, so that in the experimental setup considered here 
atoms do not easily thermalize. An alternative interesting quantity to study is 
hence given by entropy, whose value can be extracted from the density matrix 
of the system. 

Let us first discuss the fidelity of the gate in the case where the density matrix 
of the relative motion is given by the canonical thermal distribution 

Prcl (T) =2sinh-^ £ e- E ^ T \k){k\ , (18) 

where T is the temperature, ks is the Boltzmann constant, and are the os- 
cillator energies. The density operator satisfies the normalization Tr[p re i(T)] = 
1. 



7 Actually, the square sinus potential enters the single atom Hamiltonians. In this 
case it is not even possible to exactly separate the centre of mass and relative 
motion Hamiltonians. Such decoupling is indeed a peculiar property of the harmonic 
potential. However, the quartic correction can be used as a crude test of anharmonic 
effects. 
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Fig. 4. Temperature dependence of the fidelity (\<p\ = tt, At\ = 10 Th , A = 0). 
Curves from top to bottom correspond to 6 = 0.05Th o , xq = 20 (solid line), 
= 0.05T ho , x = 40 (dashed line), 6 = 0.1T ho , x = 20 (dotted line), 6 = 0.1T ho , 
xo = 40 (dot-dashed line). 

By substituting p = Pre\{T) into Eq. (15), the fidelity of the finite temperature 
motional state is found to be 




where is given in Eq. (17). The resulting behaviour for reasonable values 
of the (purely harmonic, A = 0) trap parameters is reported in Fig. 4. Notice 
that current experimental techniques allow to cool the system much below the 
highest temperature considered in the figure (which is three times the tem- 
perature defined by the oscillator energy). In the numerical calculation, for 
each considered temperature the density matrix is constructed by including 
in Eq. (18) all the eigenstates up to at least E k = It is worth mention- 

ing that the presence of anharmonicities in the trapping potential can worsen 
the large temperature behaviour of the curves shown in Fig. 4. Indeed, higher 
oscillator states, having a larger spatial extension, are more sensitive to de- 
viations from the quadratic confinement. However, for a quartic term with 
A = 1CT 2 (see previous subsection) and the parameters of Fig. 4 this becomes 
important only for /c B T > 3 (corresponding to the mean occupation number 
k ~ 2.5). 

As anticipated above, it also useful to look at the entropy increase caused by 
the gate evolution. In fact, even if the very first initialization was perfectly 
done, i.e., p = |0)(0|, after a certain number of operations one expects the sys- 
tem to be in a superposition of excited motional states. As shown by Eq. (16), 
already after a single 2-qubit gate the system is left in a non-trivial, entangled 
state of the internal and external degrees of freedom, which will spoil to some 
extent a successive application of the phase gate. Clearly, for the gate to make 
sense this must be a small effect, which is guaranteed by fidelity values close 



16 



1 

0.5 

1 

^0.5 
= 

0.1 


0.1 




.(A). 



-(B) 



icy 



(D) 



10 

N 



15 



20 



Fig. 5. Entropy S (in units of ks) as a function of the number of cycles N for 
\4>\ = 7r and Aii = 10~ 4 Ii lo . The considered trapping potential is the purely har- 
monic (A = 0) and the remaining parameters are (A) 9 = 0.05Th o , x$ = 20, (B) 
9 = 0.05T ho , x = 40, (C) 9 = 0.1T ho , x = 20, (D) 9 = 0.1T ho , x = 40. 

to 1. Qualitatively, it is also interesting to see which is the fidelity for the 
finite entropy states arising from repeated executions of the gate operation. 
Indeed, although in practice computation would not involve always the same 
two qubits, this gives an idea of the propagation of the gate error. 

Concretely, if we start with the system in the product state |^) = \x) <8> |0), 
with the same notation as in the previous subsection, after N gate operations 
|#) is transformed into the iV-step state \V(Nt)) = c gg \gg)®Y,k a ok (Nr)\k)e-^ k+1 ^ NT + 
( c ge\ge) + c eg\eg) + c ee |ee)) (g> |0)e~ l7Vr / 2 . For simplicity, we consider the case 
c gg = 1 and construct a relative motion density matrix with the iV-step coef- 
ficients ak(Nr) 

Prci(N)=J2\aok(Nr)\ 2 \k)(k\ . (20) 

k 

The entropy is S — — kBTr\p re \(N) lnp re i(7V)], where the trace can be calcu- 
lated summing over all the eigenvalues, already explicit in the diagonal form 
of Eq. (20). 

The entropy behaviour as a function of the number of gate cycles is plotted 
for reasonable parameters in Fig. 5. In the numerical calculation the entropy 
is extracted by projecting the time-evolved state onto the harmonic oscillator 
eigenstates \n). All the eigenstates up to n = 25 are included, much more 
than required for the considered parameters, ensuring that the normalization 
is always preserved. The maximum (numerical) entropy is S/ks = Inn ~ 
3.2. Note the periodic behaviour of the entropy: its value does not increase 
arbitrarily. Indeed, after a given number of cycles one has a revival of the initial 
state (i.e., only the relative motion ground state is occupied). This is true only 
for reasonable parameters: for parameters which yield a bad gate fidelity one 
does not have this revival and the entropy increases monotonically. 

As an estimate of the fidelity due to a non-zero entropy motional state we use 
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x = 40 


F =0.9964 
F(7V = 4) =0.9956 
S(N = 4)=0.39 


F =0.9998 
F(JV = 2)=0.9998 
S(N = 4)=0.01 



Table 1 

Fidelity of a single gate execution at different entropy values for the same parameters 
as in Fig. 5. Fo is the fidelity calculated for the matrix |x)(x| ® |0)(0| (zero entropy), 
while F(N) for the matrix |x)(xl ® Prei(N). The entropy S'(iV) is the same shown 
in Fig. 5. 

Eq. (15) with p = ® Prei(N). The calculation now gives 

F = ^|a fc(iVr)| 2 F fc , (21) 

and the corresponding results for some relevant values of iV are reported in 
Table 1, showing that the fidelity is not affected very much by the considered 
process. Indeed, for the used parameters the coefficients \a 0k (NT)\ 2 entering 
Eq. (20) decay very rapidly by increasing k. In other terms, although not 
directly comparable to the canonical distribution (18), the population of the 
first few states shows a behaviour corresponding to k^T 1 (in harmonic 
oscillator units). As a consequence, we also notice that the small anharmonic- 
ity (A = 10~ 2 ) previously considered does not alter significantly the results 
presented in Table 1 for the purely harmonic case. 



5 Conclusions and outlook 

In this paper we aimed at applying the method of coherent control of two- 
body dynamics, already proposed in [11] as a means for obtaining fast, high- 
fidelity quantum gates with ions, to the case of Rydberg-excited atoms. We 
found that the method works well in a realistic situation with a simple three- 
pulse sequence that allows for obtaining a two-qubit control-phase gate with 
a fidelity bigger than 99.9% in a fraction of the trap period for presently 
achievable experimental parameters. While the ideal scheme at the basis of 
the gate is completely independent of the initial motional state of the atoms, 
the approximations employed to treat the dipole-dipole interaction introduce 
deviations which worsen the gate performance for motional states different 
from the ground state. We then analyze the effect of errors in the preparation 
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of the motional state in terms of temperature and entropy. For parameters rea- 
sonably close to the ones optimized for high fidelity, the temperature of the 
initial state is found to significantly affect the gate quality only for relatively 
large values, comparable with the oscillator energy of the trapping potential. 
In addition, the entropy of the evolved state obtained from repeated applica- 
tions of the gate to the same two qubits does not increase indefinitely with 
the number of operations, but it rather oscillates in a quite regular fashion, 
exhibiting full ground-state revivals after a few gate iterations. This is likely 
to depend on symmetries in the underlying evolution, evidenced by the trivial 
iteration of the gate on the same two qubits. Even though the latter analysis 
does not convey much information on a real quantum algorithm, which typ- 
ically involves many gates between different qubits, it would be interesting 
to study the possibility that similar recurrences emerge over more complex 
gate sequences. In this case, indeed, optimization could be done not only at 
the level of a single quantum gate, but also on a broader time scale involving 
"blocks" of subsequent gates, in the sense of tailoring each gate's parameters 
in order to minimize the system entropy after a given gate sequence. A deeper 
understanding of this aspect is still missing, as it goes beyond the scope of the 
present work, and shall be the subject of future investigations. 

Let us finally comment on a possible alternative implementation of the phase- 
space displacements which allow to acquire the geometrical phase. Indeed, 
the momentum shifts obtained from the dipole kicks could in principle be 
substituted by spatial shifts, straightforwardly realizable by displacing the 
centre of the harmonic oscillator potential, provided that one manages to 
properly design the spin dependency of the trapping potential [18]. 
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